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! Abstract 

The decoherence mechanism of a single atom inside a high-Q cavity is studied, and the results 
are compared with experimental observations performed by M. Brune et al. [Phys. Rev. Lett. 76, 
1800 (1996)]. Collision dephasing and cavity leakage are considered as the major sources giving 
| rise to decoherence effect. In particular, we show that the experimental data can be fitted very 

■ well by assuming suitable values of collision Stark shifts and dark count rate in the detector 
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A two-level atom interacting with a single cavity mode, the Jaynes-Cummings model 
(JCM) is possibly the most well studied system in quantum optics and has the po- 
tentiality to constitute the basic building block of quantum computers |2j. Among all the 
intriguing phenomena related to the JCM, the oscillation in the atomic inversion probability 
- the Rabi oscillation — plays a prominent part in quantum optics. From a fundamental 
point of view, the existence of discrete Rabi frequencies provides a direct evidence of elec- 
tromagnetic field quantization [lj. In particular, when there is a dispersion in the photon 
number, the beating of different Rabi frequencies gives rise to collapses and revivals in the 
inversion probability of the atom j3]. For example, if an atom initially prepared in its ex- 
cited state |e) evolves under the influence of a cavity field, characterized by a photon number 
distribution p n and a Rabi frequency Q, the probability of finding the atom in the ground 
state \g) at a later time t is given by 



Pe 9 {t) = -X E Pn [l - COS (2fit V / ^Tl)] , (1) 



^ n =0 

which clearly demonstrates the mentioned phenomenon. 

The first observation of quantum Rabi oscillation was made by Rempe et al. some 
years ago However, their experiment failed to obtain a conclusive result due to the 
limitation on the observation time. More recently, M. Brune et al. successfully carried out 
an experiment to observe the Rabi oscillation [5], which for the first time provided direct 
and clear evidence of field quantization inside a high-Q cavity. In their experiment, Rydberg 
atoms independently interact with a photon mode in a superconducting microwave cavity 
and undergo transitions between two atomic states with principal quantum numbers 51 and 
50 respectively. The Q-factor of the cavity mode is 7 x 10 7 , which corresponds to a photon 
lifetime of 220 fis. The Rabi frequency at the center of the cavity is Sl — 507r kHz, which 
is sufficiently fast to make Rabi oscillations observable within the cavity leakage time. The 
atoms are initially prepared in the excited state and, in addition to a background 0.8 K 
thermal field with a mean photon number n t h — 0.06, the cavity field is maintained in 
coherent states with mean photon number n varying from zero to a few photons. The 
experimental data of P eg obtained from Ref. Q], for n = 0, 0.4, 0.85 and 1.77 are reproduced 
here as the boxes in Figs. 00 Despite that the data clearly revealed the Rabi oscillations 
in P eg {t), the evolution does not conform to that predicted by Eq. (JTJ) (solid lines in Fig. [TJ). 
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Instead, it was shown that a best fit to the data is given by 



OO 

Pegit) = 7, E Pn [l - exp(-rt) COS ^V^Tl)] . (2) 
^ n=0 

The oscillations were damped in such a way that P eg {t) approached 0.5 in the long term. 
This behavior is not attributable to cavity leakage, which would instead lead to P eg {t) = 1 
asymptotically. Furthermore, the photon lifetime was actually longer than the observation 
time and the effect of cavity damping due to leakage should only play a marginal role (see 
the dot-dashed lines in Fig. [TJ. It was conjectured that decoherence effect due to collisions 
with background gas might have contributed to these damped oscillations j^. 

The aim of this report is to theoretically account for the experimental results by studying 
the effects of collision dephasing and cavity leakage. We will show in the following discussion 
that these two independent mechanisms are the main culprits leading to the discrepancy 
between the experimental data and the theoretical results given by Eq. (fTj). 

In the rotating wave approximation, the hamiltonian of a two-level atom interacting with 
an ideal cavity mode is given by (in units of h = 1): 

H = -fT-Sz + u) c a)a + (of S- + aS + ^j , (3) 

where a) and a are respectively the creation and annihilation operators of the cavity field, S z 
and S± are the pseudo-spin operators of the atomic levels, and u a and u c are the atomic tran- 
sition and cavity mode frequencies respectively. Hereafter, we will assume exact resonance 
condition such that u a = uj c . It is worthwhile to note that the vacuum Rabi frequency Q 
in the hamiltonian can be smaller than the maximum Rabi frequency fio because maximum 
coupling is attainable only if the two- level atoms are exactly located at the cavity center . 
Instead, Q will be considered as a free parameter to fit the leading few Rabi oscillations. 

To include the effects of collision dephasing jmd cavity leakage on the system, we consider 
the master equation for the density matrix p 

dp 
dt 

where 



SB: 

-i[H ,p] + L c p + L f p, (4) 



Lfp = —{2apa) — a^ap — pa* a) , (5) 
with k being the cavity leakage rate (k = 4.55kHz in the experiment), and 

L c p = 27(2S,pS, - S 2 zP - pS 2 z ) . (6) 



The operator Lt results from finite leakage at the boundaries of the cavity, whereas L c 

n □ 

describes the effect of collision dephasing characterized with an energy shift 27 p, |7J. In 
general, there is an additional term due to spontaneous decay: 

r — r r + r 

L sP = --^{S z ,p}- L ^±±Lp, (7) 

with T e and T g being the spontaneous decay rates of the upper and lower energy levels 
respectively. However, in the experiment the lifetimes of the two levels are about 30 ms 
j^j], which are relatively long compared with the atom-cavity interaction time, and this term 
is ignored hereafter. Besides, we also assume that the thermal field outside the cavity is 
negligible. In terms of the density matrix, P eg {t) is given by 

p eg(t) = J2( n '9\p\n,g) ■ (8) 

n 

In the current situation, the leakage of the cavity is small. It is legitimate to take L/p as 
the small term and Eq. (J3| can then be solved by method of perturbation. First, the density 
matrix is written as 

P = P (0) + , (9) 

where p^ is a small correction to the zeroth order density matrix p(°\ It is then readily 
shown that 

a Jo) 

H ,pW}+Lj°\ (10) 



dt 
and 

dpW 



dp (o) 

■ —1 



—1 



H ,pW}+L c pW + L f p^. (11) 



dt 

Secondly, we solve Eq. (jTOjl for the initial conditions: 

^,e|/9 (0) |n,e) = p n , 

(n,g\p^\n,g^ = (n, e\p^\n + 1, g\ = (n + 1, g\p^\n, e^> = . 

When the field is prepared in a coherent state, the photon distribution is Poissonian with 
p n given by: 

_~f) n 

Pn = e- nV -. (12) 
n! 

The explicit solution to Eq. (fTTH) is: 

0,g\p {0) \0,g) = 0, (13) 



n>,g\p i0) \n,g 



n, e\p(°'\n, e 



:Pn-i 



1 — exp(— jt) 



cos(A n _it - n _i) 



1 + exp(— 7t) 



n, e\p^\n + 1, g 



COS (p n -l 

cos(A n t - 4> n ) 

COS n 

sin A„i 



-p n exp(-jt) 
2 cos 



where 



and 



\ n = dA(n+l)9? 



r 



tan</>„ = — 



2 
A, 



Eq. 



can hence be expressed explicitly as 



1 oo 

Pe 9 {t) = -Y J Pr, 
Z n=0 



exp(-Tt) 



cos(A n t - 4> n ) 



COS ( 



(14) 
(15) 
(16) 

(17) 

(18) 

(19) 



From the results, we see that the mechanism of collision dephasing gives rise to remarkable 
effects. The Rabi oscillations are damped with modified frequencies given by Eq. (fTTj) ; and 
there is also a phase shift as given by Eq. (JTHj). The solid-lines in Fig. El show the zeroth order 
perturbation results for Q = 150.2 kHz and 7 = 19.3 kHz, which is observed to provide best 
fits to the four sets of experimental data. The results obtained by numerically solving the 
exact master equation (i.e. Eq. (J3J) are also shown there as dot-dashed lines. Interestingly 
enough, there is a fairly good apparent agreement between the zeroth order results and the 
experimental data. We consider this as a coincidence due to two counter-balancing effects, 
namely the effects of cavity leakage and dark counts in the detector used in the experiment 
0| , which have been neglected in the forgoing discussion. We will return to this point later. 

To obtain a better agreement between the perturbative and the numerical results, the 
first order calculation is carried out, yielding the results: 

sin(Aot — 20o 



o,g\p {1) \o,g) = -^p 



2 sin cj)Q — 2Vtt + exp(— 7^) - 



cos 0c 



(20) 



(n,g\p m \n,g) 



sin20 n 1 
K 4nA Vn ~ X + 4 ^ + 1 ~ ~ l ' Pn ~ l \ Kt 



-k exp(— jt) < (2n 2 + n + l)p, 



sin 2<j). 
^A 



n-l 



COS \ n -lt 



n-l 
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(2n - 1)(1 - v^Pn-i + [4n 2 + 3n + 1 - 2(4n 2 + 2n + 1) sin 2 n _i]p r , 



n(n + l)p n . 

- sm A„r 



4A n _i 
(2n 2 + 3n + 2)p n 



sin A n _i£ 



2A n 
(2n - l)x/np 



4Ar. 



sin(A n t - 20 n ) 



n— 1 



2A n _i 



Otcos(A n _it - n _i) 



(21) 



n, el/9^- 1 - 1 |n, e 

^_ / sin20 ra ^ _ (2n + 3) sm 20 n+1 ^ \ + 1 ^ + 3 ) pn+i _ ( 2n + \) Pn ] K t 

+k exp(— 7rJ < sm 2<p n cos A n r 

[ 4A n 

(2n + l)(Vn + 1 - 1) + 2 cos 2</>J p n - [n + 1 + (4n 2 + 10n + 7) cos 2</> n ] 



sin A„£ 



4A n 

(n+l)(n + 2)p n+1 . (2n 2 + 7n + 5)p n+1 

■ sm A n+ it H sm(A n+ it - 20 n+1 ) 



2A n 



+i 



4A n 



+i 



(2n + l)Vn + lp w 

2 An, 



fit cos(A n t — n ) 



(22) 



/n, e | p^ 1 -* | + 1, g\ 

a I i f T^ ~ Pn+l, 

o\/n + 111 



8VnTTfl 



exp(— 7*) < p n - (3ri 2 + 8n + 6)p n+i cos \ n t 



+ 



+ (^ 2 + 4n + 2)p n+ i tan n sin A n £ + 2(n + l)(3n + 5)p n+ i cos A n+ it 

ft 2 



+2(ra + 1) p n+1 tan0 n+1 sin A n+] t + 4(n + l)(2n + l)p n — t sin A^ > . 

A n 



(23) 



The first order results with the same parameters are shown in Fig. El In the time regime 
under consideration, the first order results give excellent approximation to the exact numer- 
ical solution that includes the leakage effect. However, as shown in Fig. |31 the introduction 
of leakage obviously also worsens the agreement between our theoretical results and the 
experimental data. As suggested in Ref. |5[, we propose that the effect of cavity damping 
might have been counter-balanced by dark counts in the atomic detector, which become in- 
creasingly important at long times because of low atomic fluxes. If it is assumed that within 
the time range of the experiment, the detection rate approximately goes exponentially with 
t, an extra factor e~ at should be introduced. In other words, 



\Peg CO] experiment ^ [-Peg (0] theory • 



(24) 
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In the time regime under consideration, this is roughly the same as assuming a linear time 
dependence of the dark count rate. In order to counter-balance the effect of cavity damping, 
we find that a should assume a value around 1.59 kHz. In other words, the proportion of 
dark counts should increase from at t = to around 13% at t = 90 fxs. In Fig. 01 the 
overlying solid lines and dot-dashed lines respectively show the first-order perturbation and 
numerical results with the effect of dark counts taken into account. It clearly shows an 
excellent agreement with experimental data. 

Moreover, the results of the first order perturbation allow us to write Eq. (JHJ) in the form 

oo 

Peg{t) = J2Pnfn{t), (25) 
ra=0 

where the functions f n (t) are obtained by grouping terms in Eqs. (fTSJl - lfrSJl and d2Dj)-(J2H|). 
Instead of evaluating P eg {t) from the given photon distribution function p n j5j, we can in 
fact search for initial states that best fit the data. With the set of functions f n {t), we 
perform a least-square best fit to obtain the optimal photon distribution p(n). The results 



in which the solid lines represent the 
namely n = 0, 0.4±0.02, 0.85±0.04 



of this inversion process are shown by circles in Fig. 
theoretical results for the initial states given in Ref. 
and 1.77 ± 0.15. It is observed that the results we obtained show a fairly good agreement 
with the theoretical values. Furthermore, we perform independent best-fits on p(n) to find 
the coherent states that give the best agreement. We found that the distributions are 
respectively best fit by coherent states with n = 0.098, 0.46, 1.19 and 1.95, signifying the 
agreement of the experimental and theoretical results. 

In conclusion, we have theoretically analyzed the evolution of a single atom inside a high- 
er cavity in the presence of collision dephasing and leakage mechanisms, and the results are 
compared with experimental observations performed by Brune et al. . Collision dephasing, 
leakage and dark counts are shown to be the major factors affecting the experimental results. 
Interestingly enough, cavity leakage and dark counts tend to produce opposite effects on the 
inversion probability. By assuming suitable values of collision Stark shift and dark count 
rate, we have explicitly demonstrated that the experimental data agree nicely with the 
theoretical prediction. 
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FIG. 1: P eg (t) vs time for (a) n = 0, (b) n = 0.4, (c) ra = 0.85 and (d) n = 1.77. Boxes represent 
experimental data reproduced from Ref. j|. Solid lines are the theoretical results for k = 7 = and 
Q. = Qq = 507rkHz. Dot-dashed lines are the results when finite cavity damping with k = 4.55 kHz 
is taken into account. 
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Peg(t) 




FIG. 2: P eg (t) vs time for (a) n = 0, (b) n = 0.4, (c) n = 0.85 and (d) n = 1.77. Boxes represent 
experimental data reproduced from Ref. Solid lines are the exact analytic results with collision 
dephasing 7 = 19.3 kHz, k = and 0, = 150.2 kHz. Dot-dashed lines are the exact numerical 
results with the same set of parameters except that k = 4.55 kHz. 
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FIG. 3: P eg (t) vs time for (a) n = 0, (bVn = 0.4, (c) n = 0.85 and (d) n = 1.77. Boxes represent 
experimental data reproduced from Ref. 5j . Solid lines are the results of the first order perturbation 
for a case with collision dephasing 7 = 19.3kHz, k = 4.55 kHz and Q, = 150.2 kHz. Dot-dashed 
lines are the exact numerical results. 
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FIG. 4: Same as Fig. |31 except the introduction of dark counts with a = 1.59 kHz. 
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FIG. 5: Photon statistics for the four experimental situations. Circles represent the results ob- 
tained by inverting the optimal p n from the experimental data. Dot-dashed lines are Poissonian 



distributions that best fit the circles. Solid lines are the Poissonian distributions specified in Ref. 
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